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, Hadron production at all pT in heavy-ion collisions in the framework of parton re- 

^ combination is reviewed. It is shown that the recombination of thermal and shower 

^ ' partons dominates the hadron spectra in the intermediate px region. In d -|- Au col- 

c ■ 

lisions, the physics of particle production at any rj is basically the same as at = 0. 

• 1-^ - 

X 

' The Cronin effect is described as a result of the final-state instead of the initial-state 

interaction. The suppression of Rcp at high r] is due to the reduction of the soft parton 
density on the deuteron side, thus resulting in less pions produced by recombination, 
an explanation that requires no new physics. In Au + Au collisions large p/n ratio is 
obtained because the thermal partons can contribute to the formation of proton more 
than they do to the pion. 
PACS numbers: 23.75.-q, 25.75.Dw 
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The conventional approach to hadronization at high pt is by use of the fragmentation 
model, which has been highly successful in leptonic and hadronic collision processes. How- 
ever, the application of such a hadronization scheme to heavy-ion collisions has encountered 
a number of difficulties. There are at least five areas where the predictions disagree badly 
with data. The problems can all be resolved if fragmentation is replaced by recombination. 
The reason is simple. When a hard parton at large px is in the environment of coUinear soft 
partons that are abundant in heavy-ion collisions, then the recombination of a semi-hard 
shower parton in the parton jet with a soft parton cannot be ignored. Those processes turn 
out to dominate in the intermediate px region and contribute to the anomalies that cannot 
be understood in the fragmentation picture. Even at very high pt where fragmentation is 
valid, such a process can be interpreted as the recombination of shower partons in the same 
jet. Thus it is possible to frame all hadronization processes in terms of parton recombination 
— hence, the title of this talk. 

The five areas that can be listed at this point as puzzles in the fragmentation picture are 
the following. 

1 . The proton-to-pion ratio in central Au + Au collisions is greater than 1 at 2 < < 4 
GeV/c. 

2. In d-\- Au collisions Rcp for proton is greater than that for pion in a pt region even 
wider than that above. 

3. Azimuthal anisotropy as measured by V2 is greater for baryons than for mesons for 
Pt>2 GeV/c. 

4. The structure of jets produced in Au + Au collisions is different from that in p + p 
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collisions. 

5. Forward production in d+Au collisions is more suppressed in central than in peripheral 
collisions at nearly all px, contrary to the naive interpretation of the Cronin effect. 

Space does not permit adequate discussion of all of these problems here, although they 
were all summarized in the talk presented. We describe in this written report only items 1, 2 
and 5. A summary of the others can be found in Ref. [1], which contains the basic references 
on each topic. 

Starting with the latest phenomenon in item 5, BRAHMS data from d + Au collisions 
at RHIC at 200 GeV indicate that the central-to-peripheral ratio, Rcp, in the 1 < pr < 3 
GeV/c region decreases monotonically from a value ~ 1.8 at pseudorapidity 77 ~ —2 to a 
value ~ 0.5 at 77 ~ 3.2 [2]. This has led to the interpretation of a change of the physics 
responsible for the phenomena from the gold side 77 < to the deuteron side (77 > 0) [3]. 
For 7; < the enhancement of the particle yield that has been referred to as the Cronin 
effect is generally regarded as the result of multiple scattering in initial-state interaction. 
For 77 > saturation physics has been considered to be dominant, especially at large 77, 
so that there is suppression, instead of enhancement, in particle production [4]. Neither 
explanation takes into account any details about hadronization in the final state. The use of 
fragmentation model is inappropriate because of the known failure to explain the p/vr ratio. 
We discuss below how the Cronin effect can be understood in terms of parton recombination 
without any multiple scattering in the initial state [5] . Then we present an extension of that 
consideration to 77 > and show that the suppression in forward production can be well 
reproduced without the explicit introduction of any new physics [6]. 
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Following the formalism developed in Refs. [5, 7] for parton recombination at pr > 0.5 
GeV/c, we have for the inclusive distribution of pion in a ID description 

dN^ f dpidp2 , . , . , . 

/ Fqq,{pi,p2,7])R^{pi,P2,p), (1) 

dpdr] J Pi P2 

where the recombination function for forming a pion at p is R.j^{pi,p2,p) = {pi,P2/p)S{pi + 
P2 — p)- For p in the transverse plane so that px — p, the distribution dN^^/ (Ppdr], averaged 
over all 0, is 

dN^ 1 /"P , ^ , , , , 

^-psi dpiF,APi.P-PuV). (2) 

The joint parton distribution Fgg' has three components 

Fgg, = TT + TS + SS, (3) 

where T stands for soft parton distribution and S for shower parton distribution. At low pr 

the observed pion distribution is exponential, which suggests the form 

dN'^ 

nPi, V) = Pi^^ = Cif3, rj)p, exp 77)] , (4) 

where C{(3,r]) and T{(3,r]) are to be determined phenomenologically. Here, we use /3 to 
denote centrality, e. g., /5 = 0.1 for 0-20% centrality. We do not rely on any model to derive 
the properties of the soft partons, but determine T{(3, rf) from the data on soft pions. Our 
prediction is on the behavior in the intermediate pt region where hard partons can cause 
measurable effect. 

The distribution 5 is a convolution of the hard parton distribution rf) with transverse 
momentum k and the shower parton distribution (SPD) Si {z) from hard parton i to semi- 
hard parton j 

Sj(pi,v)-E[ dkkMk,rj)Si(p,/k) , (5) 



where kmin is set at 3 GeV/c, below which the pQCD derivation of fi{k,r]) is invahd. A 
power-law parametrization of fi{k, rj) for d + Au collisions is given in Ref. [6], and of Si {z) 
in Ref. [8]. 

At midrapidity, 77 = 0, the particle production rate at a fixed high in d + Au collisions 
increases with centrality, a phenomenon that has been known for nearly thirty years [9] , and 
referred to as the Cronin effect. The traditional interpretation of the centrality dependence 
is in terms of multiple scattering of projectile partons by the target nucleus before the 
production of a minijet by a hard scattering, which is followed by the fragmentation of 
the hard-scattered parton. Since the fragmentation process occurs outside the nucleus, the 
ratio of the production rates of proton to pion should be just the ratio of the fragmentation 
functions for the two types of hadrons h, when all else are kept the same. It is known 
that is much smaller than . However, the early data from PHENIX showed that BFcp 
for proton is greater than R^p for pion in the range 1 < pr < 3 GeV/c [10], where R^p is 
the ratio of central to peripheral production of hadron h. That finding puts into question the 
reliability of the hadronization scheme in terms of fragmentation. The fragmentation model 
has been successful in describing hadron production in leptonic and hadronic collisions. But 
in nuclear collisions its inadequacy points to the relevance of the soft partons associated 
with the nuclear medium that are absent in the collision of simpler systems. The dynamical 
process of hadronization that involves the soft partons is recombination. 

The recombination formalism summarized in Eqs. (l)-(5) can be applied to the produc- 
tion of mesons. For proton production one only needs to generalize the 2-parton distribution 
Fqqi to 3-quark distribution Fund with a corresponding generalization of the recombination 
function R to take into account the wave function of the proton in terms of the valons [5, 11]. 
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In Fig. 1 we show the results of our calculations of Rcp for pion and proton; they agree very 
well with the data [10]. The physical reason for the success is that when soft partons are 
considered in the formation of hadrons in the final state, their abundance enhances the for- 
mation of proton more than pion. Shower partons are needed to increase the hadronic pt, 
but the soft partons increase the yield. Since no initial-state broadening of the parton trans- 
verse momenta has been put in, it is the final state rather than the initial state interaction 
that is mainly responsible for the Cronin effect. Item 2 listed in the introduction as one of 
the puzzles associated with fragmentation is therefore resolved. 

The above is concerned with d + Au collisions at midrapidity. Let us now consider what 
happens at forward and backward rapidities. We can proceed to the calculation of the pion 
distribution given by (2), using (3)-(5), provided that we specify C{P,rj) and T(P,rj) in 
(4). In [5] we have determined C{P, 0) and T{j3, 0); now we need to extend them to r] > 0. 
Since the observed rapidity density at dNch/drj is an integral over the pt distribution that 
is dominated by the soft contribution at low pr, i- e., the TT term in (3), it should be 
proportional to C^(/5,?7). We can therefore determine C[j3,rj) by rescaling from C(/3,0) 



C{(5,r^)^C{(5,Q) 



(6) 



_dN,h/d7]\r,=o{P) 

Using PHOBOS data on dNch/drj{P), we can obtain C{P,rj), as shown in [6]. For T{P,rj) 
the assumption that it is independent of r] is not a bad approximation. However, a slight 
decrease with increasing 77 yields a better fit. We use 

T{f3,ri) = To{l-ePri) (7) 
where Tq = 0.208 GeV and e = 0.0205. 



With the soft parton distribution specified, we can now calculate the pion distribution for 
all pt and rj. The results for tt"*" production are shown in Fig. 2 for two extreme centralities 
in d + Au collisions and for four values of 77. Note how the distributions are progressively 
more suppressed at high pt as 77 is increased. That has to do with the hard partons at high 
k approaching the kinematical limit as rj becomes large. In addition to that suppression 
at high Pt and large r], there is also the suppression at high (3. That is a result of the (3 
dependence prescribed by (6), which corresponds physically to the fact that there are less 
soft partons in more peripheral collisions. 

A quantitative measure of the (5 and 77 dependences is the ratio Rqp given by 

„ , dNjpTdpTdr,{(5)/ {N,.m{0)) 

^''^^''^^ dN^/pTdpTdr^{(5,)/{N^u{^)) ^ > 

where the reference centrality is /3p = 0.7. We have calculated that ratio for all /? and rj 
values of the data from BRAHMS [2] . The results are shown in Fig. 3 and compared with 
the data. Evidently, the theoretical curves agree well with data. No new physics such as 
gluon saturation has been considered explicitly, although the data on dNch/drj used in (6) 
may in turn be described by saturation physics. For us here that part is a phenomenological 
input. The curves in Fig. 3 also take into account the efi:ect of momentum degradation that 
is responsible for baryon stopping, but it is a minor effect unworthy of extended discussion 
in this short summary here. 

Finally, let us summarize the situation with Au + Au coUisions. RHIC experiments have 
convincingly proven that the away-side jets are attenuated due to energy loss of partons 
traversing a thick hot medium, suggesting that the near-side jets are created mainly by hard 
collisions close to the near-side surface. Since our hadronization formalism does not trace 



the space-time properties of the coUiding system and its subsequent evolution, we use a 
phenomenological parameter ^ to denote the average fraction of hard partons that emerge 
from the bulk medium to hadronize. While such a factor ^ is 1 in d + Au collision due to 
the absence of energy loss in a cold medium, it must be included as a multiplicative factor 
on the right-hand side of Eq. (5), when applied to Au + Au collisions. All other formulas are 
the same in appearance, although the functions C{P,rj),T{P,rj) and fi{k,rj) are all different 
when the colliding nuclei are changed from d + Au to Au + Au. Again, C(/9, rj) and T(/?, rj) 
are determined from low-pr pion production data, and fi{k,ri) is obtained by calculation. 
When all three contributions from Eq. (3) are considered, ^ is the only parameter that 
must be adjusted to fit the normalization of the pion spectrum at high pr- The shape is 
predicted, since the properties of the shower partons are already fixed by other considerations 
independent of the Au + Au collisions [8] . 

Fig. 4 shows the pion distribution in px at midrapidity in Au + Au at 200 GeV [7] . The 
values of ^ used to achieve the fit is ^ = 0.07. The overall shape agrees well with data 
[12]. Note that the thermal-shower recombination is dominant in the region 3 < pr < 9 
GeV/c. The shower-shower recombination in one jet is equivalent to fragmentation, and 
is important only for px > 9 GeV/c. The thermal partons therefore have a previously 
unsuspected important effect on hadron production in the intermediate pt region. The 
significance of that effect is made convincingly clear when the p/ir ratio is studied. 

As in the d + Au collision case, the proton distribution in pt can likewise be calculated 
for Au + Au collisions. When the result is used to compute the p/n ratio, we find that the 
ratio exceeds 1 at pt ~ 3 GeV/c [7]. That is shown in Fig. 5. Such a high observed ratio 
[13] has been an anomaly in the fragmentation picture. Hadronization by recombination 
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resolves the anomaly. Similar result has also been obtained by other groups using different 
implementations of the recombination model [14, 15]. 

The success of the recombination mechanism for hadronization in the intermediate pt 
region can be traced to the following simple reasons. In recombination the hadron momentum 
is the sum of the momenta of the partons that recombine. Thus the parton momenta are 
lower than the hadron momentum; the probability of finding those lower momenta partons is 
therefore higher. That is a far more efiicient mechanism than fragmentation, which requires 
a hard parton at a much higher momentum, for which a penalty is paid in yield, and then 
the application of a fragmentation function to produce a hadron at a momentum fraction 
causes another suppression in the yield. 

To conclude we have successfully described hadron production at intermediate pr in 
heavy-ion collisions on the basis of parton recombination. For d + Au collisions we have 
shown that parton recombination is sufficient to account for all the experimental features 
found in both the mid and forward rapidity regions, namely, the enhancement of 77 ~ 
and suppression at 77 > 0. As is increased from negative to positive values, no change of 
physics has been built in, i.e., from transverse momentum broadening to gluon saturation, 
both in the initial state. We have also indicated how the pion and proton spectra in Au-\-Au 
collisions can be obtained and can explain the large p/tt ratio that agrees well with data. 
Even at very high pt where the fragmentation picture is valid, hadron production can still 
be interpreted as the result of the recombination of shower partons. Since hadronization of 
partons by recombination is a process in the final stage of evolution of the partons, it is not 
in confiict with any dynamical model that correctly describes the beginning and subsequent 
evolution of those parton. Indeed, given any predicted parton spectra in any such model, 
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the work presented here provides the necessary hnk to the observed hadronic data. 

The work reported here is the result of collaboration mainly with C. B. Yang and also 
with R. J. Pries. It was supported, in part, by the U. S. Department of Energy under Grant 
No. DE-FG03-96ER40972. 

References 

[1] R. C. Hwa, in the Proceedings for the International Symposium on Multiparticle Dy- 
namics, Sonoma County, CA, July 2004, Acta Phys. Polon. B36, xxx (2005); nucl- 
th/0410038. 

[2] I. Arsene et al, (BRAHMS Collaboration), nucl-ex/0403005. 
[3] L. McLerran, hep-ph/0402137. 

[4] D. Kharzeev, E. Levin, and M. Nardi, Nucl. Phys. A 730, 448 (2004); D. Kharzeev, Y. 
Korchegov, and K. Tuchin, Phys. Rev. D68, 094013 (2003); hep-ph/0405045. 

[5] R. C. Hwa and C. B. Yang, Phys. Rev. Lett. 93, 082302 (2004); R. C. Hwa, and C. B. 
Yang, Phys. Rev. C 70, 037901 (2004). 

[6] R. C. Hwa, C. B. Yang and R. J. Fries, nucl-th/0410111, Phys. Rev. C (to be published) 

[7] R. C. Hwa, and C. B. Yang, Phys. Rev. C 70, 024905 (2004). 

[8] R. C. Hwa, and C. B. Yang, Phys. Rev. C 70, 024904 (2004). 

[9] J. W. Cronin et al, Phys. Rev. D 11, 3105 (1975). 



10 



[10] F. Matathias (for the PHENIX collaboration), J. Phys. G: Nucl. Part. Phys. 30, S1113 
(2004). 

[11] R. C. Hwa and C. B. Yang, Phys. Rev. C 66, 205204 (2002). 

[12] S. S. Adler, PHENIX Collaboration, Phys. Rev. Lett. 91, 072301 (2003). 

[13] S. S. Adler, PHENIX Collaboration, Phys. Rev. C 69, 034909 (2004). 

[14] V. Greco, C. M. Ko, and P. Levai, Phys. Rev. Lett. 90, 202302 (2003); Phys. Rev. C 
68, 034904 (2003). 

[15] R. J. Pries, B. Miiller, C. Nonaka and S. A. Bass, Phys. Rev. Lett. 90, 202303 (2003); 
Phys. Rev. C 68, 044902 (2003). 

Figure Captions 

Fig. 1. Comparison of calculated ratios for Rcp for n and p with data from [10]. 

Fig. 2. Transverse momentum distributions of tt"*" produced in d+Au collisions at different 
pseudorapidities for two centrality cuts. 

Fig. 3. Rcp for 0-20%/60-80% (filled circles and sohd hues) and 30-50%/60-80% (open 
circles and dashed lines) for four pseudorapidities. Data are from [2]; lines are the 
results of calculations in the recombination model. 

Fig. 4. Transverse momentum distribution of 7r° in central Au-Au collisions at -\/sjvAr — 200 
GeV. Data are from [12]. 
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Fig. 5. Comparison of calculated p/tt ratio with data from [13]. The difference between the 
solid and dashed lines is described in [7]. 
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